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Explosive particle production due to parametric resonance is a crucial feature of reheating in in- 
flationary cosmology. Coherent oscillations of the inflaton field lead to a periodically varying mass 
in the evolution equation of matter and gr avitational fluctu ations and often induce a parametric 



resonance instability. In a previous paper ( iep-ph/9709273) it was shown that homogeneous (i.e 



space-independent) noise leads to an increase of the generalized Floquet exponent for all modes, 
at least if the noise is temporally uncorrelated. Here we extend the results to the physically more 
realistic case of spatially inhomogeneous noise. We demonstrate - modulo some mathematical fine 
points which are addressed in a companion paper - that the Floquet exponent is a non-decreasing 
function of the amplitude of the noise. We provide numerical evidence for an even stronger state- 
ment, namely that in the presence of inhomogeneous noise, the Floquet exponent of each mode is 
larger than the maximal Floquet exponent of the system in the absence of noise. 
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I. INTRODUCTION 



Explosive particle production at the end of the period of inflation is a crucial element of inflationary cosmology. It 
determines the effectiveness of energy transfer between the inflaton (the scalar field driving the exponential expansion 
of the Universe during inflation) and ordinary matter. As was first pointed out in B and discussed in detail in |^ |^] 
and in many other papers (see e.g. @), a coherently oscillating scalar field (such as the inflaton during the period 
of reheating) induces a parametric resonance instability in the mode equations of bosonic matter fields to which it 
couples. It has recently been realized that this resonance may also amplify gravitational waves [Q, scalar gravitational 
fluctuations p|-|l(| , and - within the limits imposed by the exclusion principle - fermionic modes [Tl]JT^ | . 

Resonance instabilities are, in general, quite sensitive to the presence or absence of noise. In a real physical system 
one expects some noise in the evolution of the scalar field driving the resonance. In the case of inflation, such noise 
could be due to quantum fluctuations (the same fluctuations which on larger length scales may develop into the seeds 
for galaxies and galaxy clusters) or to thermal excitations. In both cases, the amplitude of the noise is expected to 
be small. Nevertheless, it is important to study the effects of noise on the parametric resonance instability. 

In a previous paper ]T^ ] we studied the effects of spatially homogeneous noise on the energy transfer between 
the inflaton cf> and a second bosonic matter field x- We were able to prove that - provided the noise is temporally 
uncorrelated on time scales on the order of the period of oscillation of the inflaton - the presence of noise increases the 
rate of energy transfer between tfi and x ( see a ls° 01 ^ or similar results). More specifically, the exponential growth 
rate of a particular Fourier mode \k of X is given by the generalized Floquet exponent fi^. Then, if the noise is given 
by a space-independent function q(t), the following result holds 

> Mfc(O) Vfc, (1) 

if q(t) =/= 0. In particular, this implies that the stability bands (the bands of k values for which /i(k) = 0) disappear 
completely. 

In this paper, we extend the above result to the physically more realistic case of spatially inhomogeneous noise. 
This extension involves introducing nonlinearities into the physical system. Inhomogeneous noise couples different 
Fourier modes of \- Hence, the mathematical analysis of the dynamics is much more complicated since the partial 
differential equation describing the evolution of x can no longer be reduced via Fourier transformation to a system 
of decoupled ordinary differential equations. However, by introducing an ultraviolet cutoff and by considering the 

f 



system in a finite volume, we can reduce the problem to a system of coupled ordinary differential equations. We are 
then able to show that the methods used in extend to this case. Thus, for a given cutoff theory, the result (|l|) 
can be shown to hold. When the cutoffs are removed, we obtain the slightly weaker result 

H(q) > max fe /x fe (0) = fi max , (2) 

if q(t, x) 7^ 0, where the Lyapunov exponent fi(q) is the maximal Lyapunov exponent of the inhomogeneous system. 

By Fourier transforming the configuration x( x i t) & t each time, it is possible to define a generalized Lyapunov 
exponent fJ-k(o) for each Fourier mode. Mode mixing induced by the inhomogeneous noise has the effect of distributing 
the energy which in the absence of noise would be dumped only into the original instability bands to all of the modes 
of x- Hence, it is not unreasonable to expect that a much stronger result than (Q) holds, namely that ^k{q) becomes 
larger than the maximal Floquet index [i max of the system in the absence of noise. Indeed, we find numerical evidence 
which supports the following conjecture: 

Atfc(g) > Umax Vfc, (3) 

for a general inhomogeneous noise function q(t,x). 

The outline of this paper is as follows. In the following section we introduce the dynamical system to be studied 
and review the approach of jl3| to proving ([!]). In Section 3, the mathematically rigorous generalization of the proof 
to inhomogeneous noise is outlined, relegating some technical details to a companion paper [ p^[ . This yields a proof 
of (||). In Section 4 we indicate how the main result can be seen via a perturbative approach, and in Section 5 we 
present numerical evidence for the much stronger statement of (||). 

Our results apply for both systems with narrow and broad-band resonance. For simplicity we neglect the expansion 
of the Universe; however, we do not believe that including effects of expansion would change our main conclusions. 

II. SYSTEM WITH HOMOGENEOUS NOISE 

The model we consider consists of two scalar fields 4> an d X with an interaction Lagrangian of the form 

d ~ e<j> X 2 , (4) 

where e < 1 is a coefficient which parameterizes the strength of the interaction. In our application to inflationary 
cosmology, 4> represents the inflaton field, the scalar field which drives inflation. When inflation ends, the inflaton 



begins to oscillate coherently about its ground state. Superimposed on this homogeneous oscillation is quantum or 
thermal noise q(t,x), i.e. 

4>(t,x) = p(cjt)+q(t,x), (5) 

for some periodic function p(u)t) . 

We consider <fi(t, x) to be a background field. Our goal is to study the effects of the resonant excitation of \ during 
the period when 4> is oscillating coherently (this period is called the "reheating" or more accurately the "preheating" 
period 0). The back-reaction of the excitation of x on the evolution of <fi can play an important role [0-Bl- However, 
since this is not the subject of our study, we shall neglect this effect. The equation of motion for x is 

X - V 2 X + K + frM + «(*. *))] X = , (6) 

where p(y) is a function with period 2tt. Both the function p(y) and the noise term in (^|) derive from the interaction 
Lagrangian ([|) (and are rescaled compared to (g)). 

In the absence of noise, we can diagonalize the partial differential equation (^|) by Fourier transformation. The 
modes Xk satisfy the following second order ordinary differential equation: 

Xk+ [Wfc+p(w*)]x* =0. (7) 

where ui 2 = m 2 + k 2 . This is a Mathieu or Hill equation for each k, and as is well known, there are (as a function of 
k) stability and instability bands. In the instability bands, Xk grows exponentially at a characteristic rate given by 
the generalized Floquet index If pjuj 2 <C 1, the instability bands are narrow, whereas in the case pjuj 1 > 1 wc 
have broad-band resonance. Our analysis will apply to both cases. 

In the case of homogeneous noise which was analyzed in the Fourier modes remain decoupled: 

Xk + [4 + p(wi) + ?(<)] Xk = . (8) 

It proves convenient to write the noise in terms of a dimensionless function n(t) with characteristic amplitude 1: 

q(t) = guo 2 n{t) , (9) 

where g is a dimensionless coupling constant which is proportional to the coupling constant in (Q). 

The first step in the analysis of jl3| was to rewrite Equation (||) for fixed k (the index k will be dropped in the rest 
of this section) as a first order matrix differential equation for the fundamental solution (or transfer) matrix $ 9 (r, 0) 



•«*">- (10) 

consisting of two independent solutions Xi(i;<z) an d X2(i;?) of the second order equation (||). The matrix equation 
reads 

= M(q(t),t)$ q , (11) 
with initial conditions $ 9 (0,0) = /. Here, M(q(t),t) is the matrix 

-(-M + K.) + ( K«))i)- (12) 



For vanishing noise, i.e. g(t) = 0, the content of Floquet theory is that the solution of ( |ll| ) can be written in the 
form 

$o(t,Q) = fo(i)e Ct , (13) 

where Po(i) is a periodic matrix function with period T — 2it /u>, and C is a constant matrix whose spectrum in a 
resonance region is spec(C) = {±/z(0)} with fi(0) > 0. 

The noise is supposed to describe quantum or thermal fluctuations. For the derivation which follows, it is sufficient 
to make certain statistical assumptions about q(t). We assume that the noise is drawn from some sample space f2 
(which for homogeneous noise can be taken to be = C($R)), and that the noise is ergodic, i.e. the time average of 
the noise is equal to the expectation value of the noise over the sample space. In this case, it can be shown that the 
generalized Floquet exponent of the solutions of ([ll]) is well defined by the limit 

^ = J im lo S W n U^T, (j - l)T)\\ , (14) 

where || • || denotes some matrix norm (the dependence on the specific norm drops out in the large N limit). Further- 
more, the growth rate /j,(q) is continuous in q in an appropriate topology on Q (see Appendix of [ p"3"[ for details). 

To obtain more quantitative information about [i(q) it is necessary to make further assumptions about the noise. 
We assume: 

(i) The noise q(t) is uncorrelated in time on scales larger than T, that is, {q(t; k) : jT < t < (j + 1)T} is independent 
of {q(t; k) : IT < t < (I + 1)T} for integers I ^ j, and is identically distributed., for all realizations k of the 
noise. 



(ii) Restricting the noise q(t; k) to the time interval < t < T, the samples {q(t; k) : < t < T} within the support 
of the probability measure fill a neighborhood, in C(0, T), of the origin. 

Hypothesis (i) implies that the noise is ergodic, and therefore the generalized Floquet exponent is well defined. The 
main result which was proved in is that /i(q) is strictly larger than ^t(O); 

//(<?)> M (0), (15) 

which demonstrates that the presence of noise leads to a strict increase in the rate of particle production. The 
proof was based on an application of Furstenberg's theorem, which concerns the Lyapunov exponent of products of 
independent identically distributed random matrices {^Sfj : j = 1, TV}. 

Consider a probability distribution dA on the matrices ^ € SL(2n, R). Let Ga be the smallest subgroup of SL(2n, R) 
containing the support of dA. 

Theorem 1: (Furstenberg, JT^/J Suppose that Ga is not compact, and that the action of Ga on the set of lines in 
R 2n has no invariant measure. Then for almost all independent random sequences {<f Q SL(2n,R) with common 
distribution dA, 

lim ilog(||nf =1 ^||) = A>0 . 

AT— *oo iv 

Furthermore, for given ^1,^2 £ R 2n ? then 

lim llog(( Vl ,nf =1 * jW2 )) = A 

for almost every sequence {^j}?^. 

In order to apply Furstenberg's theorem to obtain (|l5|), we start by factoring out from the transfer matrix 0) 
the contribution due to the evolution without noise; 

$ q (t,0) = $ (i,0)* 9 (t,0) = P o (t)e ct ^ q (t,0). (16) 

The reduced transfer matrix ^f q (t,0) satisfies the following equation 

% = S(t; K^ q = %\t) ( _ ?( ° k) ° \ $ (t)* 9 , (17) 

which can be written as a matrix integral equation 
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¥,(t, 0) = 1+ [ dr S(t; k)%(t, 0) . (18) 
Jo 

By properties (i) of decorrelation of the noise, the quantities ^ q (jT,(J — 1)T) are independent and identically 
distributed for different integers j, and we can apply the Furstenberg theorem to the following decomposition of 
*,(#T,0); 

* g (NT,0) = IL? =1 9 q (jT, (j - 1)T) . (19) 

We may choose for instance the vector v\ in Theorem 1 to be an eigenvector of $>o(T, 0)* = (Po(T)e )*, the transpose 
of the transfer matrix of the system without noise, with eigenvalue e^°^ T . Then the second statement of Theorem 1 
becomes 

J^logl^.^^OWl = ^log\(v 1 ,P Q (NT)e CNT * q (NT,0)v 2 )\ 

= J L ; log(e^ NT \(v 1 ,y q (NT,0)v 2 )\) (20) 

= Mo) + ^io g |( Vl ,nf =1 vi/ 9 (jT,(j-i)T)« 2 )|. 

Taking the limit N — > oo and applying the first and second statements of Theorem 1 we obtain 

fi(q) = mo) + A > MO) (21) 

which proves the main result (JlJ). Note in particular that ( [Df ) implies that the stability bands of the system without 
noise disappear when spatially homogeneous noise is added. 

III. EXACT RESULTS 

In this section we will generalize the results of [[l3| to the case of inhomogeneous noise, modulo some technical 



points which will be addressed in a separate paper |15|. To keep the notation simple, we will take space to be one- 
dimensional. However, the analysis works in any spatial dimension. The starting point is the second order differential 
equation (||) which we can rewrite as a first order matrix operator equation 

$ 9 = M(q(t,-),t)* q (22) 

with initial condition ^(O) = / for the fundamental solution matrix with kernel <fr q (t;x,y) given by 

$ (fx v) = ( X&> X >V) X2(t;x,y)\ , 
4> q (t,x,y) - \ ±l{t . Xiy) U t;x,y) J ' ^ 



consisting of two independent kernels Xi(*i x : V) an d X^(t', x, y) of the linear operator equation (g). Here, M(q(t, x), t) 
is the matrix operator 

M (^)>') = (-(-V* + mllp(t) + q (t, X )) I) ■ ( 24 ) 

The operators are taken to act on the Hilbert space H = H 1 (R 3 ) x L 2 (R 3 ), where H 1 stands for the Sobolev space 
of L 2 functions with L 2 gradients. The first factor refers to the coordinate function, the second to the velocities. 
Physically, our choice of function space corresponds to considering finite energy configurations. 

As in the previous section, we introduce rotating coordinates in order to factor out the time evolution in the absence 
of noise: 

$,(t) = *,(t)S (t), (25) 

where $o(0 is the fundamental solution matrix operator in the absence of noise, and the product implies operator 
composition. The initial conditions are ^ q (0) = /. In Fourier space this is a generalization of the matrix given in the 
previous section (|l3|). The kernel of the matrix operator ^ q (t) satisfies the reduced differential equation 

®q(t;x,y) = J S(t;x,z 1 ;q)^f q (t;z 1 ,y)dz 1 = J J $o X (i; x, z 2 ) ( ^ ® ^ $o(*! z 2> z\ Zi, y) dz x dz 2 . (26) 

The corresponding operator equation can also be written as an integral equation 

*,(t) = 1 + / dt'S(t';q)%(t'). (27) 
Jo 

We will focus on the reduced transfer matrix operator 

* V (NT) = nf =1 V q (jT, (j - 1)T) (28) 

which describes the evolution from t = to t = NT, where T is the period of the oscillating source p(i), factored into 
matrix operators ^ q {jT, (J — 1)T) which describe the evolution in each individual period. By differentiating ( p7j ) with 
respect to the noise function q we obtain the action of the transfer matrix operator in the Lie algebra of infinitesimal 
displacements r(t, x) 

8* q (T) ■ r = £ dt %\t) ( _^ >} ° ) $ (t) . (29) 
This construction will be used in the derivation of our main results. 
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The first relevant mathematical result is that the growth rate (or generalized Floquet or Lyapunov exponent) of 
the solutions of ( |22| ) is well defined by the limit 

fi(q) = lim -±- log ||nf =1 $(jT, (j - 1)T)|| , (30) 

N— >00 JVJ 

where || • |j denotes the operator matrix norm 

||*|| = su PveHlxL 4$(v)\\ H , (31) 

where the subscript H indicates that the function norm in the Hilbert space H is used. Furthermore, the growth 
rate n(q) is continuous in q in an appropriate topology on the noise sample space. Both of these results were already 
demonstrated in the Appendix of Jll| to which the reader is referred for details. 

The Lyapunov exponent of the wave equation without noise, i.e. of the wave operator *o(£) is given by 

fi = lim jlog||* (t)|| = max keR3 fi k (0) , (32) 

t — too % 

where /ifc(0) is the Lyapunov exponent of the k'th mode of the system without noise. 

We assume that the noise function q(t, x) is given by a stochastic process which is temporally uncorrelated over 
times t larger than the period T. In this case, the quantities ^ q (jT : (j — 1)T) are identically distributed uncorrelated 
random matrix operators. The main theorem of this paper is the following: 

Theorem 2: When a spatially inhomogeneous noise q(t,x;oj) is introduced into the problem, then the Lyapunov 
exponent can only increase, i.e. 

> Mo 

for almost all realization lu of the noise. 

The main idea of the proof is to reduce the infinite dimensional operator problem to a finite dimensional matrix 
problem to which Furstenberg's Theorem in the form stated in Section 2 can be applied. First, we introduce an 
infrared cutoff by putting the system in a box. This corresponds to a discretization of the Fourier modes. The 
separation of k values is denoted by Ak. Next, we impose an ultraviolet cutoff in momentum space by eliminating all 
modes with frequency larger than k max . This reduces the problem to a finite dimensional case to which Furstenberg's 
Theorem applies. We will prove that the hypotheses of the theorem are indeed satisfied. 



8 



In order to be able to take the limit when the ultraviolet cutoff goes to infinity and the infrared cutoff to zero, 
it is important to prove that the transfer matrix ^ q in the continuum system is a compact operator perturbation 
of the transfer matrix operator without noise. This implies that the perturbation can be well approximated by the 
corresponding transfer matrix operators of the finite dimensional approximate problem obtained after imposing the 
cutoffs. Since a product of compact operators is compact, it is sufficient to prove that both $o(T, a;) — <f>o (T, x) \ p —q 
and ty q (T, x) — I are compact independently (the first term describes the difference in the transfer matrix when the 
periodic part of the perturbation, namely pit), is turned on, the second term expresses the effect of non- vanishing 
q(x,t)). 

Theorem 3: The operators $>o(T, x) — &o(T, x)\ p= q and ^ q (T, x) — I are compact. 

The compactness of $q(T, x) — &o(T, x)\ p —q follows via the known asymptotics of the eigenfunctions of Sturm- 
Liouville operators. The compactness of ^f q (T,x) — I can be shown explicitly [ p^[ . 

Given the infrared and ultraviolet cutoffs introduced above, it is natural to work in Fourier space. In this basis, 
the partial differential equation (|^) for x(t, x) becomes a system of In coupled ordinary differential equations for the 
Fourier components Xfc(t) and their momenta Xfc(t), where n is the number of independent Fourier modes. In other 
words, the transfer matrix equation ( [22] ) becomes a In dimensional vector equation with a 2n x 2n transfer matrix 
&q(t). If we order the basis vectors such that 



Xki 

Xk n 
\ Xk„ J 



(33) 



where fci, k n label the finite set of momenta we are considering, then ^o{t, x), the transfer matrix in the absence 
of noise, is block diagonal 

/*o,fe \ 

*o(t,a:) - ... . (34) 

V $ ,fc„ / 

Here, the matrices ^o.fci are the 2x2 matrices for the fcj'th mode in the absence of noise (see Section 2). The matrix 
S of (]2^), on the other hand, mixes the blocks. If we define a new matrix Q by 

S(t,x;q) = Sofea;)- 1 ^*, a:) (35) 

and denote the Fourier transform of q(t, x) by q(t, k), then 
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Q(t,x;q) 

where Qij is the 2x2 matrix 

Qij (j>i fci 

We are now ready to sketch the proof of the main theorem. It is sufficient to show that the hypotheses of the 
Furstenberg Theorem (see Section 2) are satisfied. Let us denote 

vf, = * q (jT, (j - 1)T) G SL(2n, R) . (38) 

and consider Ga, the smallest subgroup of SL(2n, R) containing the support of cL4, the set of transfer matrices which 
can be obtained from realizations of the noise drawn from the noise sample space. We must show that 

(i) Ga is not compact. 

(ii) Ga acting on the set of lines in R 2n has no invariant measure, i.e. there is no subspace of R 2n left invariant 

under the action of all of the transfer matrices corresponding to noise drawn from the sample space. 

To show that Ga is not compact, it is sufficient to consider the subset Ha C R 2n corresponding to homogeneous 
noise and to show that Ha is not compact. Thus, we take 

q(t,ki-kj) = q(t)S{ki - kj) . (39) 

Thus, Ha is block diagonal and does not mix the basis vectors corresponding to different values of k. Since the set 
{q(t)} is not compact, it follows that within each block, Ha is non-compact. 

Turning to the second hypothesis of Furstenberg's Theorem, it follows from our previous study of homogeneous 
noise (see Appendix of Jl3| ) that no subspace of R 2n corresponding to a fixed value of k is invariant under the action 
of Ha and hence of Ga- To show that there cannot be an invariant subspace which is different from a subspace of 
an R 2 associated with an individual value of k, we argue by contradiction. It is convenient to do the analysis at the 
level of the Lie algebra associated with the group of transfer matrices, namely the set of infinitesimal displacements 
given by (p9[), and then to use the implicit function theorem to recover the corresponding result at the group level. 

Assuming that there exists an invariant subspace X C R 2n , we construct a special noise function q(t, x) such that 
the associated transfer matrix takes some vector v G X into a vector w 2" X. For example, let us assume that the 
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/ Qll Ql2 ■■ Qlr, 
Q21 Q22 ■■ Q2r, 

\ Qnl Qn2 ■■ Q III 



" kj) = 











q(t, ki - kj) 



(36) 



(37) 



subspace X contains a vector v € associated with the a'th Fourier mode, but that it is orthogonal to a vector 
w G R 2 ^ associated with the /3'th Fourier mode. Consider now a noise function which is defined by 

q(t, kt - kj) = q(t)8(ki - k a )S(k 3 - hp) • (40) 



Under this noise function, the vector v is mapped onto a vector v' which is not orthogonal to w. The generalization 
of this argument to arbitrary vectors v and w will be given in ]lq | . 

IV. BORN APPROXIMATION 

In this section we follow our previous work Jl3| closely. The inflaton field is taken to be a real scalar field and is 
denoted by <j). In the period immediately after inflation, (j)(t) is assumed to be oscillating coherently about (one of) its 
ground state(s). As in the previous section, we shall, however, include a small aperiodic, inhomogeneous perturbation 
(cf. eq. (§)). 

We shall (following |l|-|3| ) take 4> to be coupled to a second scalar field x which represents matter via an interaction 
Lagrangian which is quadratic in \, for example of the form (|4|). To simplify the discussion, we neglect nonlinearities 
in the equation of motion for \. Such nonlinearities may be important and lead to an early termination of parametric 
resonance. This topic has recently been discussed extensively, see e.g. [^,^|, but is not the topic of our paper. Instead, 
we are interested whether the presence of noise such as included in Eq. will affect the onset of resonance. 

For simplicity, we will neglect the expansion of the Universe. In jl] || it has been shown that it is possible to include 
the expansion without any problems and that it does not affect the results concerning the parametric resonance 
instability. Since the equation of motion for x( x i t) is linear, we can write down the evolution equation for the Fourier 
modes of \i denoted by Xk(t) or, equivalently, by x(k, f). The equation of motion for the \ field is given by (|^), 
where (?(i,x) represents the noise which we consider as a perturbation of the driving function p(t). Taking the Fourier 
Transform of equation (1q) and using the convolution theorem we obtain the equation of evolution for the fc-mode 



Xk 



[<4 +p(u>t)] Xk = -J d 3 kq(k' - MMk'.t) (41) 



where ui\ = mi + k 2 , and Xk = x(k, t) is the Fourier transform of x(x, t). 

Note that all modes on the right hand side of Eq. (^d]) will influence the time evolution of the particular fc-mode 
we are considering. This is a very important aspect of our analysis for the case of inhomogeneous noise. 
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A straightforward calculation shows that ([IT]) is equivalent to an integral equation (resembling Volterra's equation), 
namely 

Xk(t) = x[ h) + X ( k 1] = X ( k h) - f dt'G(t, f ) J d 3 k'q(k> - k, t) X (k' 7 1) (42) 

where U is the time at the beginning of reheating and ls ^ ne general solution of the "homogeneous equation" 

X + [k 2 + m\ + p(ut)] x = 0. (43) 

For the particular case when g(k' — k, i) = q(t)5(k' — k) we obtain the equation discussed in 
Since piuit) is a periodic function, xi can ^ e written in the Floquet form 

X { k h) = e^piit) + e-^pzit) (44) 

where pi(t) and p 2 (t) are periodic functions (which include the arbitrary constants determined by initial conditions) 
and [ik is the Floquet exponent. 

The solutions XiW = e fMkt pi(t) and X2 — e~ fikt P2(t) are linearly independent and we can choose them in such a 
way that their Wronskian W — Xi(t)X2(t) — Xi{t)X2{t) is time independent (see e.g. Jlj]]). Then, the kernel G(t,t') 
in Eq.(|4^) is given by 

(45) 

Our method, in this section, consists of finding an approximate solution of the integral equation ([42]), by reducing 
it as closely as possible to the equation in the case of homogeneous noise previously obtained in In the first 

approximation we replace x(k',i) in the momentum-space integral in ( (42| ) by the homogeneous solution x^ (k', t). 
The integral then reduces to 

J d 3 k'q(k' - k, t)pi(¥, t)d Lkt + J d 3 k'q(k' - k, t)p 2 (k', t)e-^ kt (46) 

If we restrict to the case of narrow resonance, the main contribution to the above integral comes from the first 
resonance band, since higher bands are of the order o(e 2 ) or higher. In this case the integrals in the above expression 
can be replaced by integrals over a finite domain, namely, the first resonance band which, in momentum space, is 
given by the shell k\ < \k'\ < For our analysis we do not need to know the values of fci and k 2 - The expression 
( ff(f ) above reduces to 
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|k|=fc 2 f\k\=k2 

d 3 k / <z(k'-k,*)pi(k',*)e« i + / d 3 k^(k' - k,t)p 2 (k' ^^e'^ (47) 

If we assume that the periodic function p2 has an amplitude of the same order as (or less than) the function p\ 1 
then we can neglect the integral containing p2, since its integrand is exponentially decreasing. Thus we are left with 
only the first integral in (|47|). Now, for this integral, we apply the following Mean Value Theorem p8[ : 

Theorem (Mean Value in 3 Dim): Let f,g:Vc 1Z 3 — ► 1Z be continuous functions, where V is a compact domain 
in 1Z 3 . Then there exists &x £ V such that 



/ /(x) 5 (x)d 3 x = /(x) f ff (x)d 3 x. 
Jv Jv 



Using this theorem, ((47|) reduces to 



k|=fc 2 



k|=jfei 



d 3 k'q(k' - k,i>i(k',i) = e^A(k, t) , 



(48) 



(49) 



where k is inside the resonance band. This is our final expression which estimates the contribution of all modes in 
the evolution equation of the mode k. So, using the above result, our earlier evolution equation ( pi~i| ) reads 

Xk + K + (p(wt)] Xk = ~A(k, t)e^ (50) 

A closer analysis of the function A{k,t), defined in ((49|), shows that it is bounded. This is because p\ is periodic 
and the noise is by assumption a random function with a small amplitude which is identically distributed over time 
periods of T . Thus, we can use the approximation 

\A{k,t)\ < M^ik^e, (51) 

where Mi and au are, respectively, upper bounds of the functions pi(k',i) and q(k' — k, t) as both variables k' and 
t are varied. The last two factors in ( |5l|) represent an estimate of the volume of the first (and dominant) instability 
band. Also, let be the upper bound on P2(k', t). 



Wk c ) 



H(k 4 ) 
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FIG. 1. The Floquet exponent fj, is plotted (vertical axis) against k/uj (horizontal axis), for a general case of the Hill 
equation 

In the resonance band [ki, k?], we must consider the Floquet exponent as a function of k. The exact formula for [i^ 
can be complicated, but it can be approximated very well by a circle of arc centered at the midpoint of the interval 
[ftijfe]- Let k c denote the central mode with the maximum exponent. Note that fJ-(ki) = /^(ta) = (see Fig. 1). 
Then, if k ^ k c , we have \x-^ < /ife, for all k satisfying \k — k c \ < \k — k c \. We will now consider the subinterval 
[&3) fei] C [ki, fc 2 ] "centered" at k c for which this is the case, and we shall show that for values of k in this subinterval 
the noise does not destroy the parametric resonance. 

From Eqs. (fl|) and (|§) it follows that 

X? (t) = - f dt'G(t, t')A(k, t')e^' (52) 



where the Green's kernel G(t, t') is given by (45). Using ( |5l| ) it can easily be shown that 

<a k f(k,At)e^ At (53) 

where At = t — ti and the function f(k, At) is given by 

/(k, At) = M?M 2 W- 1 e ttkU (l - e'^ kAt )euj 3 — (54) 

Comparing j5^ ) with the growth rate of the homogeneous solution xj^, that is e MfeAt , we find 

^M<^mAt) (55) 

for all k in the subinterval fc 3 < |k| < k^ inside the resonance band (a shell in the 3-dim momentum space) [fci, ^2]- 

The above inequality shows that small amplitude noise (which can be controlled by the parameter <7fc) will not 
change the exponential rate of growth in the parametric resonance regime. 

This method applies only under very restrictive hypotheses such as the assumption of narrow resonance band 
dynamics, applicability of the First Born Approximation method, and yields only a result in a subinterval of the 
resonance band which may be much narrower than the band itself (in fact, for k — k c , only for the central mode k c we 
could guarantee that noise does not eliminate the exponential growth). However, based on the numerical results to 
be presented in the following section we gain some confidence that the methods could be generalized to less restrictive 
hypotheses (perhaps including the case of broad resonance). This will be the subject of further studies. A very 
important aspect in our analysis is the fact that the noise considered here is quite general. This is good since there 
are many possible sources of noise during and before the end of the inflationary period. 
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V. NUMERICAL RESULTS 



In order to estimate the effects of inhomogeneous noise on the resonant field through numerical approximations we 
take the "discrete" Fourier Transform of equation @ to get 

Xk(t) + [u 2 k +p(u>t)] xk(t) = -5^&'-k(*)Xk'(*) » (56) 

k' 

where x( x ^) = Ek Xk(t) exp[ik.x], etc. 

Let us analyze the time evolution of a given mode k. We can write ([56|) in the following form 

Xk(<) + K + J>(wt) + &(*)] Xk(t) = ^k , (57) 

where i^k = — Ek'^k ?k'-k(i)Xk'(*) does not contain Xk and then may be thought of as an external force. 
The general solution to <p7\j has the form 



Xk = x£ + Xl, (58) 
where Xk is the solution to the homogeneous equation 

Xk (t) + H + p{ut) + q Q (t)] Xk(t) = , (59) 
and Xk 1S the particular solution of ( |57j ) 

Xi{t) = f G{t,t')F^t')dt' . (60) 
G(t,t') is the Green function associated to ( |59| ) 

^ xi(*')x2(*)-xi(*)x2(*') , . 

G{t,t) = ^7 , (61) 



where xi and X2 are the two independent solutions of (£>9|) and W is the Wronskian, which in the present case is a 
constant. 

Equation (J5^) is the same as the evolution equation for the k-mode of the boson field coupled to the inflaton during 
the pre-heating phase of the Universe and in the presence of homogeneous noise. This problem was studied at length 
in According to the results obtained in that analysis, the presence of ergodic noise strictly increases the rate of 
particle production. In other words, the generalized Floquet exponent in the presence of noise Hk{qo) is greater than 
/tifc(O), the Floquet exponent for q {t) = 0. We can write Hk{q ) — Mfc(0) + Xk{q ) with \k(q ) > 0. We then have two 

15 



distinct situations to analyze depending upon whether the k mode is in a resonance band or not. Modes in a stability 
band have /Ltfc(O) = 0, while the modes in a resonance band have /ife(0) > 0. 

To simplify the discussion we neglect the contribution of q (t), i.e. set Xk(q ) = in Eq. (|5^). This assumption will 
not change the qualitative results derived below. Moreover, for the situations we are interested in, where the noise is 
in general small when compared to the amplitude of the inflaton, and for modes close to the center of the resonance 
band, Xk(q ) is small compared to /ifc(0) and can be neglected. 

Let us first assume that k is in a resonance band. Then according to the theory of the Hill equation, the solution 
to (^) with q (t) = is 

Xi{t) = e"*<°>* Pl (t) + e-^p^t) (62) 

where p\ and P2 are bounded (periodic) functions of time. The solution to (^) for the modes within any resonance 
band and for late times can be written as 

Xk (t)~ e w(°)* Pl (t)+^(t). (63) 

From the above result we see that Xk(^) has an exponentially growing component independent of XkW- Therefore, the 
exponential growth of a mode in a resonance band can only be destroyed by inhomogeneous noise if there is an exact 
cancellation between XkW an d Xk(0- This is only possible if XkW * s gi ven by (at least for late times) —e^^piit), 
which is very unlikely since the initial conditions for all the modes are involved (see below). Notice that XkW can m 
fact grow exponentially in time and reinforce the increase of Xk(t), since all modes are coupled through the term on 
the right hand side of (|56|). The result is a stimulated resonance effect on the particular components, even for modes 
in a stability band (see below). 

We now study the time evolution of non-resonant modes. Our goal is to estimate the rate of growth of such modes. 
For the sake of simplicity, we take q (t) ~ as above. Therefore, the homogeneous solution XkW i s gi yen by two 
independent periodic functions Xi(t) and X2(^)- As we are interested in the exponentially growing components, we 
are left with only the particular solution XkW t° compute. In the following we proceed to find an upper bound for 
this "particular" solution. The results hold for all modes in a stability band. In order to estimate this solution, which 
is given by (|60|), we first need an approximation for fk(t). Let us then explicitly determine the contribution from the 
homogeneous and particular solutions of each mode. Using (|58|) and the definition of Fk(t), we can write 
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F k (f) = - «*-*(*) (e^WPvW + e-^^QvM+xUt)) , (64) 

where -Pk(i) and Qk(t) are periodic functions of Let us recall that /J,k(Q), the Floquet exponent, is (according to 
our definition) zero for modes in a stability band and positive for resonant modes. To begin with, we neglect bounded 
and decaying components in the above equation and keep just exponentially growing terms. Thus, using ( |6l| ) and 
recalling that in the present case the solution \i and \i that enter the Green's function G(t, t') are periodic functions 
of t, we can write 

k'GB^*" k'#k^*° 

+terms not growing exponentially in time , (65) 

where B stands for the sub-sets of the k-space R 3 , such that k is in a resonance band. Continuing our approximations, 
we consider the fastest growing term on the right had side of (|6|) . Let k c denote such a mode so we can cast ([35]) in 
the form 

X P k (t) = - f G(t,O9k c -k(Oe^ (0)t >k c (i) - J2 t GMgvMWvV) + R, (66) 

J U k'/k J *° 

where R stands for terms that grow slower than e /ifcc '- ^ t . The first term in the above equation guarantees an exponential 
growth of X\J\t) with the same rate as the fastest growing resonant mode. It is also worth noting that this leading 
term does not depend on the initial condition of the mode itself, but just on the initial conditions of the mode k c . 
This result holds for every mode in a stability band and is confirmed by numerical calculations as depicted in Figure 
2(a). 

Note that in equation ( |66|) k c =^ k, since we assumed that k is in a stability band while k c is the central mode of 
the broadest resonance band of the specific model considered. In a narrow resonance band regime k c corresponds to 
the central mode of the first resonance band. In a broad resonance regime, however, k c may not be the central mode, 
but it is still the mode which has the fastest exponential growth. For instance, in the model studied in [Q, k c = is 
the mode with the greatest Floquet exponent. 

To verify the previous qualitative results we performed a numerical analysis of the problem. Following the notation 
of Ref. |H, 

we choose p(t) = Acos(2t), where A is a constant, and the noise given in the form qk(t) — g \ n(t)m{k). 
Here, g is a positive constant, n(t) is a random function of time with characteristic time T — T^ 1 and amplitude 
equal to unity, and m(fc), also a random function with characteristic rate and amplitude 1, depends only upon the 
magnitude of k. We use the dimensionless time t ( = ojt) so that the evolution equation for \k reads 
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Xk + [El + Acos(2t)] Xk = \gn(t)^2xk>m k <- k , (67) 

k> 

where A and El are normalized with respect to the frequency u of the inflaton field. Let us stress that the amplitude 
of the noise was chosen to be very small, in such a way that all the corresponding homogeneous samples of the noise 
would not change the solution significantly. This choice assures that the effects we are about to show are consequence 
solely of the inhomogcneity of the noise. 

Figure 2 shows a particular case of broad resonance where El is of the form E\ = A + k 2 (k is also in units of ui), 
with Eg = A = 30. In this particular case, nik is chosen to be a smooth function of k which helps us to show that the 
growth of modes outside a resonance band does not depend upon their initial conditions (see below). 





gLi 1 1 1 1 1 1 1 1 1 1 

1 1.01 1.02 1.03 1.04 1.05 1.06 1.07 1.08 1.09 1.1 

(b) 

FIG. 2. The maximum amplitude of Xk(t) in t € [0, 40] is plotted (vertical axis) against E^/X (horizontal axis), E„ — A = 30. 

(a) No noise is present, g — 0. 

(b) The noise is characterized by Ft = 10 and g — 0.3 (small amplitude and large noise frequency). 

This figure shows clearly the stimulated resonance of the modes in a stability band. Fig. 2(a) shows a particular case 
of broad resonance in the absence of noise. The random peaks arise because we have chosen random initial conditions 
for both Xk(t) and Xk{t) at t = 0. In Fig 2(a) the ratio between the amplitudes of the modes in the stability band 
(for El/\ > 1.05) and the central mode is O(10~ 4 ), while in Fig 2(b), which shows the results when inhomogeneous 
noise is present, that ratio is O(10 _1 ). Observe that in Fig 2(b) all modes outside the resonance band have almost 
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the same amplitude. The absence of peaks in that region demonstrates the independence of the amplitudes on the 
initial conditions of such modes. Those amplitudes depend strongly on the initial conditions and Floquet exponent 
of the central mode Xk c {t)- Only 100 modes were taken into account in this example. 

(a) 




FIG. 3. The Floquet index fi(k) is plotted (vertical axis) against k 2 /Eq (horizontal axis), E 2 — X = 30. 

(a) No noise is present, g = 0. 

(b) The noise is characterized by Ft = 10 and g — 0.1 The noise is homogeneous (m fc /_ fe = 8(k' — k)). 

(c) Same parameters as in (b) . The noise is inhomogeneous. 

Figures 3(a)-3(c) show even more explicitly the effects of inhomogeneities in the noise. In this case k is normalized 
in such a way that it assumes 400 different values between feo = end k max = 2_Eo- In 3(a) no noise is present and we 
can identify the first resonance band. Figure 3(b) shows the result when we introduced homogeneous noise. The noise 
is random in both t and k and the amplitude is small, g = 0.1. Very small effects due to the noise are observed. Only 
modes within the resonance band grow exponentially, exactly as in the case of Fig. 3(a). On the other hand, in Fig. 
3(c) we observe a very distinct behavior of the modes outside the resonance gap. All the modes have explosive growth, 
and we cannot distinguish the stability region from the resonance band anymore. The inhomogeneity of the noise has 
large effects even in the case studied here, where the amplitude of the noise is small compared to the amplitude of 
the potential p(t). 
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FIG. 4. The same as figure 3(c), but here with Ak twice smaller than in 3(c). 

Another interesting point is related to the abovementioned spacing Afc between two successive modes. We can 
easily see how the Floquet exponent depends on Afc by comparing Fig. 3(c) with Fig. 4. In the last case we used 
the same parameters as in 3(c), only the normalization of k is different: 400 modes arc placed between fc m j n = 
and k m ax — V%Eq. rrik'-k wa s chosen to be a symmetric (otherwise random) matrix in k and k 1 in both the cases 
of Fig. 3(c) and Fig. 4. The Floquet exponent increases when the number of Fourier modes is increased. This can 
be understood since it corresponds to decreasing Afc between two neighboring modes. As a result, we are adding the 
contribution of a larger number of very close modes which are inside the resonance band, all of them contributing 
with exponentially growing amplitudes, on the right hand side of Eq. (|S7[), 

VI. CONCLUSIONS 

In this paper we have shown (modulo some technical points which will be demonstrated in a companion paper) 
that spatially inhomogeneous noise in the oscillating field which induces the parametric resonance instability does not 
decrease the Floquet exponent of the instability. This extends the results of [Q where it was shown that homogeneous 
noise leads to an increase in the Floquet exponent for each Fourier mode. Crucial to our proof was the truncation of 
the dynamical system to a finite dimensional one by means of infrared and ultraviolet cutoffs, use of the Furstenberg 
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Theorem concerning the Lyapunov exponent of products of identically distributed random matrices, and the use of 
compactness arguments to carry over the main results as the cutoffs are removed. 

We have also provided numerical evidence (backed up with approximate analytical calculations) which indicate 
that a much stronger result is true, namely that inhomogeneous noise will spread the exponential growth of the most 
resonant mode of the system without noise to all of the Fourier modes. An interesting challenge is to provide a 
mathematically rigorous proof of this stronger result. 
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